Breakdown of the Onsager reaction field theory in two dimensions 
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It is shown that for the spin 1/2 anisotropic Heisenberg model the result for the transition 
temperature is completely insensitive to the anisotropy within the Onsager reaction field theory, 
which yields a vanishing T c in two dimensions, in total conflict with the expected finite T c , increasing 
with anisotropy. This establishes that the breakdown of the Onsager reaction field theory in two 
dimensions is not limited to the Ising model, but actually extends over the whole range of anisotropy 
< \ Jij/Jij\ , \ Jfj/Jij\ < !• Therefore, for the isotropic case, the related results of vanishing T c and 
exponential dependence of spin correlation length should be seen as mere coincidence. 
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The reaction field theory first applied by Onsager to 
the many dipole problem, Q and subsequently extended 
to interacting spin systems, involves an improvement 
over the local mean field theory by subtracting from the 
orienting field a reaction field, involving that part of the 
polarization which is due to correlations. Thus, in the 
context of spin systems, the effective orienting field on a 
given spin excludes that part of the spin polarization of 
surrounding spins which is due to their correlation with 
that spin. It was also found M that "the Onsager pre- 
scription of subtracting out the reaction field is the neces- 
sary modification of molecular mean field theory in order 
to guarantee the fluctuation-response function relation of 
statistical mechanics, i.e., susceptibility = fluctuation." 
The inclusion of the reaction field introduces a degree of 
freedom which allows for self consistency between suscep- 
tibility and correlation. 

Viewed in terms of an improvement over mean field 
theory, the reduction in the magnetic transition temper- 
ature T c , relative to the mean field transition tempera- 
ture T C MF , has been discussed for the isotropic Heisenberg 
model in three dimensions. || While the original work on 
the Ising model Q does not discuss the transition tem- 
perature, formally the treatment is identical to that for 
the Heisenberg model, resulting in an identical final ex- 
pression for the correction in T c : 



J(q)/J(Q) 



(1) 



Here J(q) = J2j Jije lCi ' rij is the Fourier transform of the 
spin coupling J^- between spins at sites i and j, and Q 
is the ordering wavevector. 

The implications of this identical result for the Heisen- 
berg and the Ising model, with regard to the continuous 
spin-rotation symmetry, have not been discussed in the 
literature so far. This is especially pertinent in two di- 
mensions, where the q sum in Eq. (1) diverges logarith- 
mically for short-range spin interactions, yielding a van- 
ishing transition temperature. Hence, for the Ising model 
there is obviously a breakdown of the Onsager method in 
two dimensions, as even in the absence of a continuous 



spin-rotation symmetry, a vanishing T c is obtained, in- 
dicating no connection between the Onsager theory and 
the Mermin- Wagner theorem. 

In view of this breakdown for the Ising model, is 
the Onsager-theory result of vanishing T c for the two- 
dimensional Heisenberg model, then, in accordance with 
the Mermin- Wagner theorem ? In other words, is the 
vanishing T c result crucially dependent on the continu- 
ous spin-rotation symmetry, or is it a mere coincidence 
? This is the question we will explore in this paper. 
For this purpose, we consider an anisotropic spin 1/2 
Heisenberg model with anisotropy in spin space. With 
the removal of the continuous spin-rotation symmetry, 
if a non-vanishing T c is obtained, then indeed the On- 
sager method is sensitive to the continuous symmetry. 
However, if the vanishing T c result persists, then this 
will confirm that the continuous symmetry has no role to 
play, and the vanishing T c is a coincidence, indicating a 
breakdown of the Onsager reaction field method in two 
dimensions. 

These questions are relevant in view of the recent use of 
the Onsager method in quantitative studies of the Neel 
temperature Tn, and the spin correlation length £, for 
the spin 1 /2 Heisenberg antiferromagent in two and three 
dimensions. Q In addition to the vanishing Tn, an ex- 
ponentially large spin correlation length at low temper- 
atures was also obtained in two dimensions. This ex- 
ponential dependence arises from the same logarithmic 
divergence in Eq. (1) which leads to the vanishing Tn, 
hence the two results are not independent, but follow 
from the same source. 

We find, not surprisingly, that the introduction of the 
anisotropy in spin space does not affect the vanishing 
T c result within Onsager theory, indicating a complete 
insensitivity to the continuous symmetry. Whereas, an 
anisotropy driven gap in the spectrum of transverse spin 
fluctuation modes (magnons), and a non- vanishing T c of 
the order of the gap are characteristic features of the 
anisotropic Heisenberg model. The inescapable conclu- 
sion then is that the Onsager method breaks down in two 
dimensions, and the results of vanishing Tn and expo- 
nential spin correlation length for the Heisenberg model, 
following from Eq. (1), are coincidental. 
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I. ANISOTROPIC HEISENBERG MODEL — 
ONSAGER CORRECTION 

We consider the spin 1 /2 anisotropic Heisenberg model 
with anisotropy in spin space, 

H = — "^^[JijSixSjx + jfjSiySjy + J^SizSjz] , (2) 

m 

where = Jfj and | Jfj \ > J?- 1 , Jfj \ , so that order- 
ing first occurs, with decreasing temperature, in the z- 
direction. For an antiferromagnet the spin couplings are 
all negative. Considering the spin-symmetric, paramag- 
netic state at T > T c , in the presence of an infinitesimal 
external magnetic field in the z-direction, the Weiss mean 
field, 



H 



Weiss 



i 



(3) 



acts as an effective orienting field on the spin at site i 
within the Weiss mean-field theory, and the resulting po- 
larization (Siz) is self-consistently determined from 

(Su) = £ Si^ Si * HT ^l £ ■ ( 4 ) 

In the limit of infinitesimal field, upto linear terms 

(&*> = [£ Sl/(2S+l)].f3H^ iss 
{Su} 

/ c2 \ a tt Weiss 

— Wiz/free spin •i~'- n i z 



~ ,0 u Weiss 
X - H iz 



(5) 



where 



X O = 0(Sl) ireesp - m =pS(S+l)/3. (6) 



This result is actually valid for any spin quantum num- 
ber S, for which Si Z takes all integer-step values between 
—S and S. 

From Eqs. (5) and (3), we obtain after Fourier trans- 
formation, 

<S z (q)>=E<^).e iq - r * 



= x u [#r(q) + ./*(q){^(q)>], 



(7) 



where J*j = £^ , J z (q') exp — iq'.(rj — rj). This yields 
the RPA susceptibility as the magnetization response to 
an infinitesimal external field, 



X*(q) = (S,(q)>/flr(q) 



i-x°J»(q) 



(8) 



The mean field transition temperature T c MF is obtained 
from the divergence of the spin susceptibility at the or- 
dering wavevector Q, and from Eqs. (8) and (6), 



Tr = uQ).(sth 



cc spin 



J z (Q)S(S + l)/3, (9) 



which is valid for general spin S, and exhibits no effect of 
the anisotropy because the analysis involves an isolated 
single spin in the frozen mean field. 

Going over now to the Onsager reaction field correc- 
tion, an improvement over the mean-field approximation 
is obtained by subtracting a term from the mean field 
which is proportional to the mean field itself: 



H 



off 



£ J*i 

3 



(S jz ) 



£ 

3 



^ij J%j 



) 



(10) 



where the parameters Ay are temperature dependent, 
and characterize the Onsager reaction field. Ay is related 
to the correlation (Si Z Sj z ), and the subtracted polariza- 
tion Xij(Si Z ) is that part of the polarization (Sj Z ) which 
is due to correlation with the spin at site i. 

Following the earlier analysis for the spin susceptibility 
leading to Eq. (8), one obtains 



X 



l-X%/»(q)-A) 



(11) 



where A = Ylj^ijJfj- As before, T c is obtained from 
the divergence in the spin susceptibility at the ordering 
wavevector Q, and one obtains 



rpMF 



= [1-\{T C )/J Z {Q)] . 



(12) 



The reaction field reduces the net mean field, and hence 
the transition temperature. This shows that an impor- 
tant characteristic of the Onsager correction is that it 
acts within the mean field theory. This is in contrast 
to the spin-wave theory which involves transverse spin- 
fluctuation corrections about the mean-field state. ^| 

The presence of the A term in the susceptibility ex- 
pression in Eq. (11) introduces a degree of freedom 
which allows for self-consistency between the susceptibil- 
ity and correlation. The fluctuation-dissipation theorem 
provides the following connection between the zero-field 
spin correlations and the spin susceptibility: || 



/3<s«^> = £xs(qK q - ry 



(13) 



Now, for the isotropic model in the paramagnetic state 
at T > T c , all integer-step spin states for Si Z between —S 
and S are equally probable, and therefore for zero exter- 
nal field (Sf z )° = S(S + l)/3 for general spin S. For the 
anisotropic model higher \Si Z \ values are more likely due 
to the anisotropy, and therefore (Sf z )° = a.S(S + l)/3, 
where a > 1 is a temperature-dependent factor which 
accounts for the anisotropy. 
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However, for the spin 1/2 case a = 1, as the rela- 
tion (Sf z }° = S{S + l)/3 (with S* = 1/2) continues to 
hold even in the presence of the uniaxial anisotropy, be- 
cause both spin states Si Z — ±1/2 are symmetrical and 
therefore equally likely in the paramagnetic state. This 
also implies that (Sf z )° — (Sf z )f lccsp i n . It is this equal- 
ity which is responsible for identical final expressions for 
the Onsager-corrected T c for the Ising, Heisenberg and 
anisotropic Heisenberg models. Therefore, from Eq. (6), 



(3(S? Z )° =paS(S+l)/3 = a X c ' 
which leads to the sum rule: 



(14) 



*V*(q) - A) 



(15) 



Choosing A such that the above sum rule is satisfied en- 
sures the self-consistency between the susceptibility and 
the correlation determined from it. Taking the limit 
T -» T+ in Eq. (15), using X°(T C )J Z (Q) = T C MF /T C 
from Eqs. (6) and (9), and eliminating A from the above 
equation and from Eq. (12) determining T c , we obtain 
after straightforward algebra the following expression for 
the corrected transition temperature, 



T 

± c 



1 



^(q)M(Q) 



(16) 



For the S = 1/2 anisotropic model a = 1, as dis- 
cussed earlier, so that the result for the transition tem- 
perature is totally insensitive to the anisotropy. The (log- 
arithmic) divergence in two dimensions, and therefore 
the vanishing T c result, remain unchanged even with the 
anisotropy. This is in total conflict with the expected 
result for the anisotropic model of a finite T c , increas- 
ing with anisotropy. This establishes that the break- 
down of the Onsager reaction field theory in two dimen- 
sions is not limited to the Ising model, but also holds for 
the anisotropic Heisenberg model in the whole range of 
anisotropy < \JfJJ^\ < 1. Therefore, the vanishing 
T c result for the isotropic case should be seen as a mere 
coincidence. 

It is instructive to examine the origin of the diver- 
gence. The divergence is at the wavevector q = Q, and 
follows from the definition of the transition temperature 
where the susceptibility must diverge. Thus (not sur- 
prisingly) it is present for the Ising model, as well as for 
the anisotropic Heisenberg model, and has nothing to do 
with the presence of any soft Goldstone mode connected 
with the continuous spin-rotation symmetry. 

Continuing with the spin 1/2 case, for which the sum- 
rule simplifies to X) 9 Xz( t l) = X°i from Eqs. (13) and 
(15) one obtains: 



S(S + l)/3 ' 



(17) 



confirming that the Onsager parameters are related to 
the spin correlations, and providing for a self-consistent 
determination of the correlation function (Si Z Sj Z }° in Eq. 
(13). 

For vanishing T c , it follows from Eq. (12) that A = 
J 2 (Q), and hence |Ay| = 1, indicating perfect NN spin 
correlations at T = 0. Thus T c vanishes simply because 
the reaction field exactly cancels the mean field term in 
Eq. (10), reducing the effective orienting field to zero. 
This is a very different mechanism from that involved 
in a two-dimensional system possessing continuous spin- 
rotational symmetry, where the vanishing of T c is associ- 
ated with the long-wavelength Goldstone modes. 

In conclusion, it has been shown that for the 
anisotropic S — 1/2 Heisenberg model the Onsager reac- 
tion field theory yields a vanishing T c in two dimensions 
in the whole range of anisotropy < | / J* | , | Jf, / J* | < 
1, and generally shows a total insensitivity to the 
anisotropy factor. This result is in total conflict with 
the expected finite T Cl increasing with anisotropy, and 
confirms that the breakdown of the Onsager reaction 
field theory in two dimensions is not limited to the Ising 
model, but actually extends over the whole range of 
anisotropy. Therefore, for the isotropic case, the related 
results of vanishing T c and exponentially large spin cor- 
relation length should be seen as a mere coincidence, un- 
related to the Goldstone mode associated with the con- 
tinuous spin-rotational symmetry. 

The Onsager reaction field theory clearly underesti- 
mates T c , as seen dramatically in two dimensions. On 
the other hand, the neglect within this theory of trans- 
verse spin fluctuations must result in an overestimate of 
T c . It is interesting to note that for hypercubic lattices 
these two effects exactly compensate each other, resulting 
in another accidental coincidence. The Onsager theory 
result in Eq. (1) for the transition temperature exactly 
coincides with the Neel temperature Tn for the quantum 
Heisenberg antiferromagnet on a hypercubic lattice, as 
obtained within the Green's function theory, Q and the 
self-consistent spin fluctuation theory. || 
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